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Abstract 

Relative  Motion  in  a  system  of  coupled  rigid  bodies  can  yield  global  reorientEition  (or 
phase  shift).  We  give  a  formula  to  compute  such  a  phase  shift  and  interpret  the  same 
in  geometric  terms.  The  theory  of  connections  in  principal  bundles  provides  the  proper 
setting  for  questions  of  the  type  addressed  in  this  paper.  A  related  optimal  control  problem 
leads  to  singular  riemannian  geometry. 

1.  Background 

In  1987,  Jair  Koiller  introduced  me  to  the  work  of  M.V.  Berry  [i]  and  B.  Simon  [15] 
on  geometric  phases  in  classical  and  quantum  physics.  Inspired  by  his  remarks,  I  obtained 
formulas  for  analogous  geometric  phases  in  multibody  systems.  It  became  apparent,  from 
conversations  with  Jerrold  Marsden,  that  principal  connections  were  involved.  Around  this 
time,  A.  Shapere  was  completing  his  thesis  on  gauge  kinematics  of  deformable  bodies  under 
F.  Wilczek  [12].  In  the  hands  of  Marsden.  Montgomery  and  Ratiu,  the  gauge-theoretic 
links  between  phases  and  the  reconstruction  problem  in  hamiltonian  mechanics  [10]  [6]  [S] 
have  become  very  clear.  In  the  present  paper,  I  give  a  ’’bare  hands”  derivation  of  the  phase 
shift  formula  for  coupled  planar  rigid  bodies,  formulate  an  optimal  control  problem,  and 
present  the  geometric  picture  that  generalizes  the  example  of  this  paper. 

My  understanding  of  the  gauge- theoretic  framework  for  reconstruction  phases,  owes 
a  great  deal  to  the  many  conversations  I  have  had  on  this  subject  with  Jerrold  Marsden, 
Richard  Montgomery  and  Tudor  Ratiu.  Thanks  are  also  due  to  A1  Shapere  and  Frank 
Wilczek  for  stimulating  discussions  and  references.  This  work  is  an  off-shoot  of  the  program 
to  understand  the  dynamics  of  multibody  systems  begun  in  collaboration  with  Marsden, 
Sreenath  and  others. 

*  This  work  was  supported  in  part  by  the  AFOSR  University  Research  Initiative 
Program  under  grants  AFOSR-S7-0073  and  AFOSR-90-0105,  and  by  the  National  Science 
Foundation’s  Engineering  Research  Centers  Program:  NSFD  CDR  S3  03012  and  also  by  the 
Army  Research  Office  through  the  Mathematical  Sciences  Institute  at  Cornell  University. 


2.  Planar  N-body  Dynamics 

Consider  a  system  of  planar  rigid  laminae  connected  by  single  degree  of  freedom  pin 
joints.  The  bodies  are  free  to  float  and  move  freely  in  space  (in  the  absence  of  external 
forces  and  ignoring  self-collisions).  See  Figure  1. 


Figure  1 

The  instantaneous  configuration  of  such  a  system  relative  to  the  center-of-mass  of  the 
assembly  is  given  by  the  N-tuple  of  absolute  orientations  of  frames  rigidly  affixed  to  each 
body,  say  at  its  center-of-mass.  We  denote  it  as  q  =  (51,  •••  ,?iv)  €  TW,  the  N-torus. 

Let  u>i  A  qi  denote  the  absolute  angular  velocities  and  u>  =  (uq,  •  •  •  ,u>n)T ■  The  joints 

are  assumed  to  be  articulated  via  motors  with  drive  torques  ut,  i  =  1,  2,  *  -  - ,  (iV  —  1). 
It  can  be  shown  [16]  [18]  that  the  N-body  system  is  governed  by  a  Hamiltonian  control 
system: 

Q  =  W1  P 

P  =  ~Yq  T  0*’  J^)_1  P)  +  ^  ^ 

where  J  is  a  configuration-dependent  N  x  N  matrix  that  defines  the  free  Hamiltonian, 

H  =  \  (P,  r\q)  p),  (2.2) 

and  the  N  x  (N  —  1)  matrix  B  satisfies, 

f  1  *  =  J 

Bij  =  -1  i  =  j  +  1,  (2-3) 

l  0  otherwise. 

The  matrix  J  depends  only  on  the  differences  of  absolute  angles  and  is  hence  invariant 

under  the  5 1  action,  (</i,  •••  ,qx)  —>  (<7i  +  <?h  <72  +  <t>  » 1  •  • ,  <7Ar  +  (j>),  where  <f>  £  S1 . 

N 

The  angular  momentum  fi  =  Yj  pi  is  conserved  (control  torques  are  internal). 


The  space  of  equivalence  classes  of  configurations  under  the  diagonal  S1  action  on 
T'V  is  referred  to  as  the  (labelled)  shape  space.  It  is  the  space  T'v  1  of  joint  angles 
9i  —  qi+ 1  -  qi,  i  =  1,  2,  •  •  •  ,N  -  1.  Poisson  reduction  by  5 1  yields  the  controlled 


dynamics  on  P  —  T ^  1  x  RjV, 


Bi  = 

dH 

dH 

dpi+i 

dpi 

dH 

Pi  = 

d9 !  + 

U\ 

dH 

dH 

P2  = 

d9\ 

+  de  2 

i  =  1,  2,  •••  ,1V-  1. 

(2.4) 

+  u2 


PN- 1  — 

PN  = 


a# 


dH 

d&N-i 
UN- 1- 


UN— 2  +  UN  —  1 


For  any  choice  of  controls,  the  motions  are  confined  to  the  symplectic  leaves  p  =  J_l  Pi  = 

»=i 

constant,  in  P.  Much  of  this  is  by  now  standard  and  may  be  found  in  [16]  [18]  [n].  For 
Poisson  reduction,  see  [7]  [5].  We  pose, 

The  Reconstruction  Problem :  Given  a  motion  in  shape  space  (symplectic  leaf)  reconstruct 
the  motion  in  configuration  space  (unreduced  phase  space). 

A  partial  solution  to  the  reconstruction  problem  can  be  given  by  considering  the  splitting : 

u  =  q  =  «i  e  +  M  9  (2-5) 


where  M  =  [Mij]  is  the  N  x  (N  -  1)  matrix  defined  by, 


f  0,  i  =  1 

Mij  -  <  !,  *  >  i 

l  0,  otherwise, 


and  e  =  (1,  1,  •  •  • ,  1)T  is  a  column  vector  in  .  Then,  using  p  =  e  •  Jw,  we  get, 


O-’l 


P 

e  •  Je 


e-J  M6 
e  ■  Je 


(2.6) 


Integrating  both  sides  of  (2.6)  over  a  path  7  in  the  shape  space,  we  get  a  formula  for  the 
absolute  reorientation  or  phase  shift  of  the  body  numbered  1: 


The  first  term  in  the  formula  (2.7)  depends  on  the  time-parametrized  path  7 (<)  via  the 
centrifugal  potential  energy  and  we  refer  to  this  as  the  dynamic  phase.  The  second 
term  in  (2.7)  depends  only  on  the  path  and  not  on  its  parametrization  and  is  referred  to 
as  the  geometric  phase.  The  factor  e  •  Je  appearing  in  (2.7)  is  always  a  positive  quantity 
and  is  simply  the  locked  inertia  of  the  body  at  the  shape  9.  Thus  when  p  —  0,  only  the 
geometric  phase  remains  and  it  corresponds  to  a  retrograde  rotation  of  the  reference  body 
1. 


Example  (2  body  problem) 

Let  p  =  0.  Then, 
Aqi  =  A5f°metrlc 


=  -  /  — 
J  h  + 


2  +  £  di  d2  cos  ( 9 ) 

I2  +  2  e  d\  d2  cos  ( 9 ) 


where, 


ii  =  It  +  £  d],  i  —  1,2 


m  1  m2 


mi  +  m2 

mi  and  J,-  are  respectively  the  mass  and  inertia  about  the  body-center-of  -mass  of  the  ith 
body  and  di  are  as  in  Figure  2. 

For  interpretations  of  such  phase  shifts  in  the  context  of  gymnastics  and  diving  examples 
and  for  explicit  integration  of  (2.8),  see  the  paper  of  Frohlich  [3].  Formulas  analogous  to 
(2.7)  are  known  for  the  simple  rigid  body  in  3  dimensions.  SeefgjfG]. 

For  N  >  2,  a  closed  curve  7  can  be  the  boundary  of  a  smooth  surface  in  shape  space. 
Then,  by  Stokes’  theorem,  the  geometric  phase  is  given  by, 

^geometric  =  _  ^  d  (2.9) 

r 

where  T  is  any  surface  in  5  with  boundary  cT  =  7. 

3.  An  Optimal  Control  Problem 

Substituting  (2.6)  in  (2.5),  we  get 

r  u  1  L  e  cT  Jl  .  • 


_£ _ ' 

J  e 


e  +  1  - 


M  6, 


=  Xn(<i)  +  D(q)v, 


(3.1) 


Figure  2 

where  v  =  9,  the  shape  velocity,  is  treated  as  a  control.  Immediately,  we  can  ask: 

(a)  Is  the  system  (3.1)  accessible? 

(b)  Can  we  “solve”  the  variational  problem, 

<i 

min 

„(■) 

to 


J  <  v,  v  >s  dt 


subject  to  the  conditions, 

q(io)  =  q° 

=  31  ? 

If  q°  and  qx  are  in  the  same  equivalence  class  in  shape  space,  then,  (b)  is  just  the  isoholon- 
omy  problem  of  Montgomery.  Here  <  ,  >s  is  a  riemannian  metric  on  shape  space. 

If  fi  =  0  and  accessibility  holds,  then  one  derives  a  metric  geometry  on  configuration 
space  from  the  solution  to  (b).  This  is  a  singular/  subriemannian/  nonholonomic  geometry 
in  the  sense  of  Brockett  [2],  Hermann  [4],  Strichartz  [19],  and  Vershik  -  Gershkovich  [20]. 

Letting  fj.  =  0,  and  using  the  flat  metric  on  the  shape  space,  the  necessary  conditions 
for  an  optimal  control  are  that  there  exist  a  state-costate  pair ,  p{t))  satisfying 

q  =  D(q)  DT(q)p 
i>  =  -  ^  (1  PT  O(q)  DT(q)pj 
V  =  DT(q)p 


(3.2) 


See  Sreenath  [17]  for  a  study  of  these  equations  for  N  small. 

4.  The  Geometric  Picture 

Everything  I  have  done  above  is  part  of  a  general  geometric  picture  first  recognized 
by  Montgomery,  Wilczek  and  Shapere.  I  give  a  simple  treatment  below. 

Let  ( Q ,  K,  V,  G)  be  a  simple  mechanical  system  with  symmetry  i.e.  Q  is  the  con¬ 
figuration  space,  K  is  the  kinetic  energy /riemannian  metric,  V  :  Q  — ►  1R  is  a  G-invariant 
potential  where  G  acts  freely  on  Q  by  isometries.  The  system  may  be  extended  to  include 
controls  via  the  Lagrange  -  D’Alembert  principle. 

A  basic  object  of  interest  is  the  principal  bundle 


Q 

i 

Q/G  =  S 

where  S  is  the  shape  space.  Let  Q  denote  the  Lie  algebra  of  G  and 

(?q  :  Q  -»  TQq 

£  |->  (q  (?) 

define  the  infinitesimal  action.  Let  Ig  :  Q  xQ  —*  1R  define  a  symmetric  bilinear  form, 
(£,  v)  =  aq  l)  f°r  3-11  £■,  V  €■  Q.  We  also  define  1“  :  Q  — >■  Q*  by  setting 

ijj  =  a*  aq  and  7\b  :  TQq  — ►  T*Qq  is  the  usual  Legendre  transform.  It  is  an  easy 
exercise  to  verify  that 

J*  :  TQ  -»  Q * 

Wq  H-+  Wq) 

defined  by  setting  J *  (u>g)f  =  (K* w q)(^Q{q))  is  the  conserved  momentum  map  of  the 
free  hamiltonian  system.  Admissible  controls  are  internal  controls,  i.e.  they  also  leave  J # 
invariant.  Consider  the  splitting, 

TQq  =  ( Vert)q  ©  ( Hor)q 

Wq  -  aq  I*  1  H  +  ( Wq  -  CTq  ij  fJ.) 

where  71  =  J#(wq)  •  The  splitting  has  the  equivariance  property  w.r.t.  the  G  action 

and  defines  a  principal  connection.  It  appears  that  this  principal  connection  was  known 
to  Smale  and  Rummer.  One  can  define  a  (/-valued  1-form 


A  :  Q  ->  T*Q  ®  Q 

Aq  :  TQq  t— »  Q 


W  Ht  A/(5,  w) 

=  «a-v,)-v;a-V 


The  holonomy  of  this  connection  captures  the  geometric  phase.  In  the  example  of  this 
paper  the  splitting  (3.1)  defines  a  connection.  The  Lie  algebra  =  R  and  hence  the 
connection  form  is  just  an  R- valued  form.  Note  that  if  f  :  Q/G  — >  Q  is  a  cross  section  of 
the  bundle  Q  — *  Q/G  then  one  can  pull-down  A  to  A  =  f* A  :  Q/G  -4  T*Q  ®  Q.  In 
our  example,  since  the  bundle  over  shape  space  is  trivial,  one  has  cross-sections  and  the 
pull-down  version  of  the  connection  1-form  is  precisely  given  by  the  form 

e  •  J  M  dO 
e  •  J  e 


in  equation  (2.7). 

Wilczek  and  Shapere  refer  to  A  as  the  master  gauge  field  [14],  [13]. 

5.  Final  Remarks 

In  a  complete  manuscript  under  preparation  we  discuss  other  aspects  of  the  problem 
of  this  paper  including  the  accessibility  question  (recall  Ambrose  -  Singer  theorem)  in 
geometric  terms. 
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